HYDRODYNAMIC LIMIT FOR A BOUNDARY DRIVEN STOCHASTIC LATTICE 
GAS MODEL WITH MANY CONSERVED QUANTITIES 
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Abstract. We prove the hydrodynamic limit for a particle system in which particles may have different 
velocities. We assume that we have two infinite reservoirs of particles at the boundary: this is the so-called 
boundary driven process. The dynamics we considered consists of a weakly asymmetric simple exclusion 
process with collision among particles having different velocities. 



1. Introduction 

Interacting particle systems have been the subject of intense studies during the last 30 years due to the 
fact that, in one hand, they present many of the collective features that are found in real physical systems, 
and, in the other hand they are, up to some extent, mathematically tractable. Their study has led in many 
cases to a more detailed understanding of the microscopic mechanisms behind those collective phenomena. 
We refer to p3] for further references, and to [S] for recent results. 

Since their introduction by Spitzer [21] . the simple exclusion process and the zero-range process have 
been among the most studied interacting particles systems, and they have served as a test field for new 
mathematical and physical ideas. 

In the last years there has been considerable progress in understanding stationary non equilibrium states: 
reversible systems in contact with different reservoirs at the boundary imposing a gradient on the conserved 
quantities of the system. In these systems there is a flow of matter through the system and the dynamics 
is not reversible. The main difference with respect to equilibrium (reversible) states is the following. In 
equilibrium, the invariant measure, which determines the thermodynamic properties, is given for free by the 
Gibbs distribution specified by the Hamiltonian. On the contrary, in non equilibrium states the construction 
of the stationary state requires the solution of a dynamical problem. One of the most striking typical 
property of these systems is the presence of long-range correlations. For the symmetric simple exclusion this 
was already shown in a pioneering paper by Spohn [22) . We refer to [H [7] for two recent reviews on this 
topic. 

The hydrodynamic behavior of the one-dimensional boundary driven exclusion process was studied by 
[8], [9] and [15]. Also, Landim, Olla and Volchan [18] considered the behavior of a tagged particle in a 
one-dimensional nearest- neighbor symmetric exclusion process under the action of an external constant, and 
made connections between the behavior of a tagged particle in this situation and a process with infinite 
reservoirs. 

We consider a stationary non-equilibrium state, whose non-equilibrium is due to external fields and/or 
different chemical potentials at the boundaries, in which there is a flow of physical quantities, such as heat, 
electric charge, or chemical substances, across the system. The hydrodynamic behavior for this kind of 
processes in any dimension has been solved by [HUH]. Nevertheless, they have solved this problem only for 
the case where the unique thermodynamic observable quantity is the empirical density. 

Our goal is to extend their results to the situation when there are several thermodynamic variables: density 
and momentum. It is not always clear that a closed macroscopic dynamical description is possible. However, 
we show that the system can be described by a hydrodynamic equation: fix a macroscopic time interval 
[0,T], and consider the dynamical behavior of the empirical density and momentum over such an interval. 
The law of large numbers for the empirical density and momentum is then called hydrodynamic limit and, 
in the context of the diffusive scaling limit here considered, is given by a system of parabolic evolution 
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equations which is called hydrodynamic equation. Once the hydrodynamic limit for this model is rigorously 
established, a reasonable goal is to find an explicit connection between the thermodynamic potentials and 
the dynamical macroscopic properties like transport coefficients. The study of large deviations provides such 
a connection. The dynamical large deviation for boundary driven exclusion processes in any dimension with 
one conserved quantity has been recently proved in [llj . 

The dynamical large deviations for the model with many conserved quantities is studied at [12] , and the 
hydrodynamic limit obtained in this article is important for such large deviations. 

The model which we will study can be informally described as follows: fix a velocity v, an integer 
N > 1, and boundary densities < a v (-) < 1 and < (3 V {-) < 1; at any given time, each site of the set 
{1, . . . , N — 1} X {0, . . . , N — is either empty or occupied by one particle at velocity v. In the bulk, 

each particle attempts to jump at any of its neighbors at the same velocity, with a weakly asymmetric rate. 
To respect the exclusion rule, the particle jumps only if the target site at the same velocity v is empty; 
otherwise nothing happens. At the boundary, sites with first coordinates given by 1 or N — 1 have particles 
being created or removed in such a way that the local densities are a v (x) and (3 v (x): at rate a v (x/N) a 
particle is created at {1} x {x} if the site is empty, and at rate 1 — a v (x) the particle at {1} x {x} is removed 
if the site is occupied, and at rate (3 v (x) a particle is created at {N — 1} x {x} if the site is empty, and at 
rate 1 — (3 v (x) the particle at {N — 1} x {i} is removed if the site is occupied. Superposed to this dynamics, 
there is a collision process which exchange velocities of particles in the same site in a way that momentum 
is conserved. 

Similar models have been studied by [TJ [TUl HD]. In fact, the model we consider here is based on the 
model of Esposito et al. [TU] which was used to derive the Navier-Stokes equation. It is also noteworthy that 
the derivation of hydrodynamic limits and macroscopic fluctuation theory for a system with two conserved 
quantities have been studied in [3 . 

Under diffusive time scaling, assuming local equilibrium, it is not difficult to show that the evolution of 
the thermodynamic quantities is described by the parabolic system of equations 



where v = (1,V\, . . . , Vd), p stands for the density and p = (pi, . . . ,pd) for the momentum. F is a thermo- 
dynamical quantity determined by the ergodic properties of the dynamics. 

Therefore, the purpose of this article is to define an interacting particle system whose macroscopic density 
profile evolves according to the partial differential equation given by with initial condition 



with D being a suitable domain, and the equality on the boundary being on the trace sense. 

This equation derives from the underlying stochastic dynamics through an appropriate scaling limit in 
which the microscopic time and space coordinates are rescaled diffusively. The hydrodynamic equation (jl.ip 
thus represents the law of large numbers for the empirical density and momentum of the stochastic lattice 
gas. The convergence has to be understood in probability with respect to the law of the stochastic lattice 
gas. Finally, the initial condition for (jl.ip depends on the initial distribution of particles. Of course many 
microscopic configurations give rise to the same initial condition (pcbPoX')- 

The article is organized as follows: in Section [2] we establish the notation and state the main results of 
the article; in Section [3l we prove the hydrodynamic limit for the particle system we are interested in; the 
proof of a Replacement Lemma needed for the hydrodynamic limit is postponed to Section 21 in Section [5] 
we prove the uniqueness of weak solutions of the hydrodynamic equations also needed for the hydrodynamic 
limits. 



Let Tff = {0, . , . ,N — l} d = (Z/NZ) d , the d-dimensional discrete torus, and let D d N = S N x T d N \ with 
Sn = {1, ■ ■ ■ i N — 1}. Further, let also V C R d be a finite set of velocities v — (vi, . . . , Vd)- Assume that V 
is invariant under reflexions and permutations of the coordinates: 




(1.1) 



(p,p)(Q, ■) = (po,P o )(0 and iP,P){t,x) = (p,P)b(x),x£ dD, 



2. Notation and results 



(«i, . . . , Vi-i, — Vi, v i+1 , ...,v d ) and (uo-(i), ■ ■ ■ , ?V(<i)) 
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belong to V for all 1 < i < d, and all permutations a of {1, . . . ,d}, provided (v\,... ,Vd) belongs to V. 
Finally, denote the d-dimensional torus by T d = [0, l) d = (M/Z) d . 

On each site of Dfj, at most one particle for each velocity is allowed. We denote: the number of particles 
with velocity v at x, v G V, x G Dfj, by r)(x, v) G {0, 1}; the number of particles in each velocity v at a site 
x by r\ x = {i](x,v):v G V}; and a configuration by i] — {rj x ;x G Dfj}. The set of particle configurations is 

X N = ({0,1} V )^ . 

On the interior of the domain, the dynamics consists of two parts: (i) each particle of the system evolves 
according to a nearest neighbor weakly asymmetric random walk with exclusion among particles of the 
same velocity, and (ii) binary collision between particles of different velocities. Let p(x, v) be an irreducible 
probability transition function of finite range, and mean velocity v. 

xp(x, v) = v. 

X 

The jump law and the waiting times are chosen so that the jump rate from site x to site x + y for a particle 
with velocity v is 

1 d 1 

3 = 1 

where 5 x , y stands for the Kronecker delta, which equals one if x = y and otherwise, and {ei, . . . , e^} is the 
canonical basis in R d . 

2.1. The boundary driven exclusion process. Our main interest is to examine the stochastic lattice gas 
model given by the generator £ n which is the superposition of the boundary dynamics with the collision 
and exclusion: 

£ N = N 2 {£ b N + £%+£%}, (2.1) 

where £ b N stands for the generator which models the part of the dynamics at which a particle at the boundary 
can enter or leave the system, £ C N stands for the generator which models the collision part of the dynamics 
and lastly, £ e ^ models the exclusion part of the dynamics. Let / be a local function on X^. The generator 
of the exclusion part of the dynamics, £ e ^ , is given by 

(£^/)(»?) = E E v{x,v)[l-r,(z,v)]P N (z-x,v)[f{rf'*' v ) -f(r,)}, 

where 

i](y, v) if w = v and z = x, 
rj x ' y ' v (z,w) = ^ v( x i v ) if w = v an d z = y, 
r)(z,w) otherwise. 



We will often use the decomposition 

'-N = *~n "f" > 



nex pcx,\ . pex,1 



where 



and 



(^•V) (v) = \ E E v)[i - r,(z, v)] w) - m\ , 

vev x ,x+zeD% 

\z-x\=l 

(£ e N 2 f) (v) = ^ E E v) [i - n(z, v)] P (z x, v) [f( V ^n - m\ . 

vev x^x+zeDfj 

The generator of the collision part of the dynamics, £ C N , is given by 

veDf, qeQ 

where Q is the set of all collisions which preserve momentum: 

Q = {q = (v, w, v' , w') G V 4 : v + w = v' + w'}, 



the rate p(y, q, rf) is given by 

p(y,q,v) = v(y,v)v(y,w)[l - r)(y,v')][l ~ 7}(y,w% 

and for q = (vo,Vx,V2,V3), the configuration rj y ' q after the collision is defined as 

y,qi \ J T](y, Vj+2) if z — y and u — Vj for some < j < 3, 
^ ' ' \ r](z,u) otherwise, 

where the index of should be taken modulo 4. 

Particles of velocities v and w at the same site collide at rate one and produce two particles of velocities 
v' and w' at that site. 

Finally, the generator of the boundary part of the dynamics is given by 

(£ b Nf)(v) = E E^t^^I 1 " «)1 + 1 1 ~ ".(^j^^^I/ta 1 ^) - /(,,)] 

X\— 1 

+ E E^^/^t 1 -^^^)^^-^^))^^)]!/^^)-/^)]. 

where x = (x%, . . . , aid), 



a x ' v r](y,w) 



1 — 77(2;, 10), if w — v and y = x, 
rj(y,w), otherwise. 

and for every v £ V, a v ,/3 v £ C 2 (T d_1 ). We also assume that, for every v £ V, a v and j3 v have images 
belonging to some compact subset of (0, 1). The functions a v and f3 v , which affect the birth and death rates 
at the two boundaries, represent the densities of the reservoirs. 

Note that time has been speeded up diffusively in (|2.1[) . Let {r](t),t > 0} be the Markov process with 
generator Cn and denote by {S'f , t > 0} the semigroup associated to Cn- 

Let -D(R+, Xn) be the set of right continuous functions with left limits taking values on Xn- For a 
probability measure /i on Xn, denote by the measure on the path space D(NL+, Xn) induced by {rj{t) : 
t > 0} and the initial measure /x. Expectation with respect to P M is denoted by E M . 

2.2. Mass and momentum. For each configuration £ £ {0, 1} V , denote by 2o(£) the mass of £ and by 
lk{0' k — 1; • • • ' * ne momentum of £: 

/o(o = E^)> i*(o = £«*£(«)• 

Set J(£) := (io(£), ■ ■ ■ ,ld(0)- Assume that the set of velocities is chosen in such a way that the unique 
quantities conserved by the random walk dynamics described above are mass and momentum: ^2 xeD d I(rj x ). 
Two examples of sets of velocities satisfying these conditions can be found at [10] . 

For each chemical potential A = (Ao, . . . , Ad) £ R d+1 , denote by m\ the measure on {0, 1} V given by 

m A (0 = ;J^exp{A- /(£)}, (2.2) 

where Z(X) is a normalizing constant. Note that m\ is a product measure on {0, 1} V , i.e., that the variables 
{£(«) : v £ V} are independent under m\. 

Denote by fj,^ the product measure on Xn, with marginals given by 

V\{V ■ V(x,-) = = Tn\(0, 
for each £ in {0, 1} V and x £ D N . Note that {r)(x, v) : x £ D N , v £ V} are independent variables under /i^, 
and that the measure /j,^ is invariant for the exclusion process with periodic boundary condition. 
The expectation under /x^ of the mass and momentum are given by 

p(X) := £ M £ [J (»7x)] = J>W> 
Pfe(A) := E^n [I k (r)x)} = E u fc^( A )- 



In this formula 9 V (X) denotes the expected value of the density of particles with velocity v under my. 

exp |a + Efe=i A fc«fe| 



e v {\) := E mx Uv) 



1 + exp |a +Efe=i A fe«fe} 



Denote by (p,p)(A) := (p(A),pi(A), . . . ,Pd(A)) the map that associates the chemical potential to the 
vector of density and momentum. It is possible to prove that (p, p) is a diffcomorphism onto it C R d+1 , the 
interior of the convex envelope of {!(£),£ <E {0,1} V }. Denote by A = (A ,...,A d ) :H^M d+1 the inverse of 
(p,p). This correspondence allows one to parameterize the invariant states by the density and momentum: 
for each (p,p) in il we have a product measure v^ p — p^ p p ^ on Xn- 

2.3. Hydrodynamic limit for the boundary driven exclusion process. Let D d = [0,1] x T d_1 . Fix 
po : D d — > M + and p : D d — > R d , where p Q = (po,i, • • • ,Po,d)- We say that a sequence of probability 
measures (pn)n on Xn is associated to the density profile po and momentum profile p , if, for every 
continuous function G : D d — > R and for every 5 > 0, 



lim /i 

N^oo 



and for every 1 < k < d 



N 



n ■ 



lim p 



JV 



q 



W £ G[^)l { Vx )- f G(u)p (u)du 



> S 



> 5 



0, 



= 0. 



xeD d M 

Fix T > and let (B, || • \\b) be a Banach space. We denote by L 2 ([0, T],B) the Banach space of measurable 
functions U : [0, T] -> B for which 



ic/l 



L 2 ([0,T],B) 



/ ||C/ t |||dt < oo. 
Jo 



Moreover, we denote by if 1 (£) d ) the Sobolev space of measurable functions in L 2 (D d ) that have generalized 
derivatives in L 2 (D d ). 

For x=(x 1 ,x)e {0,1} x T d_1 , let 



d(x) 



Fix a bounded density profile po ■ D d 
function (p, p) : [0, T] x D d -> R + x " 
equations 



= T/veviPvW^rfvtx), . . .,v d f3 v (x)), if xi = 1. 

1+ , and a bounded momentum profile p : D a 



(2.3) 
A bounded 



+ air is a weak solution of the system of parabolic partial differential 
dt(p,p)+E v evvlvVx(0v(A(p,p)))} = |A(p,p), 

(2 ' 4) 

(p,p)(0,-) = (A),Po)(0 and iP,P)(t,x) = d(x),xE {0,1} x T d -\ 

if for every vector valued function H : [0, T] x D d -> R d+1 of class C 1 ' 2 ([0,T] x D d ) vanishing at the 
boundary, we have 

f H(T,u)-(p,p)(T,u)du- [ H(0,u)-(p ,p Q )(u)du 

J D d J D d 

= [ dt [ du\{p,p){t,u)-d t H{t,u) + \{p,p){t,u)- V d 2 u M(t,u) 
J o J d« { 2 1<l<d 



dS b(u) ■ d Ul H(t, u) 



f dt f 

JO J{l}xT li - 1 

- / dt f du VtJ-x(0„(A(p,p))) V Vi d Ui H{t, 

J J D* vev l<i<d 



dS a(u) ■ d Ul H(t, u) 



dS being the Lebesgue measure on T d 1 . 

We say that that the solution (p, p) has finite energy if its components belong to L 2 ([0, T], H 1 (D d )): 

pT 



and 



ds 



ds 



D'l 



\\V p(s,u)\\ du < oo 



D d 



\\Vp k (s,u)\\ du) < oo, 



for k = 1, . . . , d, where V/ represents the generalized gradient of the function /. 
In Section [5] we prove that there exists at most one weak solution of the problem 

Theorem 2.1. Let ({i n )n be a sequence of probability measures on Xjy associated to the profile (po,p Q ). 
Then, for every t > 0, for every continuous function H : D d — > R vanishing at the boundary, and for every 
5 > 0, 



lim I 



l - E h 

ad 



N d 



N 



D d 



H(u)p{t, u)du 



> S 



and for 1 < k < d 



lim Pj. 



Jfd E B (^) iMt)) - J H(u) Pk {t,u)di 



> S 



= 0, 



= 0, 



where (p,p) has finite energy and is the unique weak solution of equation 



The strategy to prove Theorem 12.11 is to use a replacement lemma, together with some estimates on 
Dirichlct forms and entropies for this boundary driven process. 



3. Hydrodynamic limit for the boundary driven process 

Fix T > 0, let A4+ be the space of finite positive measures on D d endowed with the weak topology, 
and let M. be the space of bounded variation signed measures on D d endowed with the weak topology. Let 
A4+ x A4 d be the cartesian product of these spaces endowed with the product topology, which is metrizable. 

Recall that the conserved quantities are the mass and momentum presented in subsection 12.21 For k = 
0, . . . , d, denote by 7r t ' N the empirical measure associated to the A;th conserved quantity: 

Kt' N = E I k(Vx(t))S x/N , (3.1) 
xeD d N 

where 6 U stands for the Dirac measure supported on u. We denote by < tt^' N ,H > the integral of a test 
function H with respect to an empirical measure n t ' . 

Let D([0, T] ,A4 + xA4 d ) be the set of right continuous functions with left limits taking values on A4+ xAi d . 
We consider the sequence of probability measures on D([0,T],M+ x M d ) {Qn)n that corresponds to the 
Markov process = (TTt' N , ■ ■ ■ , Trf :N ) starting from pjy. 

Let V be an open neighborhood of D d , and consider, for each v £ V, smooth functions K% : V — > (0, 1) 
in C 2 (V), for k — 0, . . . , d. We assume that each K v k has its image contained in some compact subset of 
(0, 1), that the restriction of K = J2 v ev( K 0' v i K i' ■ ■ ■ > v d K d) t° {0} x T d_1 equals the vector valued function 
a(-) defined in (|2.3p . and that the restriction of n to {1} x T d_1 equals the vector valued function &(•), also 
defined in (|2.3j) . in the sense that k(x) = d(xi,x) if x £ {0, 1} x T d_1 . 

Further, we may choose k for which there exists a constant 9 > such that: 

k(ui, u) = d(—l, u) if < Ux < 9 , 

k{u\, u) = d(l, u) if 1 — 9 < u\ < 1 , 

for all u G T d_1 . In that case, for every N large enough, vf? is reversible for the process with generator C h N 
and then (—N 2 C h N f,f) y N is positive. 
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Wc then consider v£ the product measure on Xjv with marginals given by 



v K {v ■ ■) = £} = toa (k ( k ))(O i 

where m\(-) was defined in (12. 2 [I . Note that with this choice, for iV sufficiently large, we have that if 
x £ {1} x T^T 1 , then £^ [77(0:, «)] = a v (x/N) and if a; G {N - 1} x T*r\ then E u n[t](x,v)] = (3 v {x/N). 

3.1. Entropy estimates. Let us recall some definitions. Recall that is the semigroup associated to the 
generator £ N = N 2 (C e ^ + C C N + C b N ). Denote by f t — ft the Radon-Nikodym derivative of fj, N with 
respect to i/„ . For each function / : Xn — > R, let D v n{J) be 



where 



Atf (/) = D% (/) + D c u» (/) + (/), 

^1 (/) = E E E p ^ - u ) / [vT^o - Vm] 2 



and 



D k^ = H E /pfo^V/M-vTwJ "Ufa), 

D ^(/) = E E /[a»(5/i\0(l-»7(*,«)) + (l-««(2/JV))»?(*»«)]5 



*e{i}xT£ 



E E 



[p v (x/N)(l - r](x,v)) + (1 - p v (x/N)) V (x, v)]x 



Proposition 3.1. There exists a finite constant C — C(a,/3) such that 

d t H(fi N S?\^) < -N 2 D„ N {f t ) + CN d . 



(3.2) 



Proof: Denote by C* the adjoint operator of Cm with respect to . Then, f t is the solution of the 
forward equation 

d t f t = N 2 £*J t , 



Thus, 



d t H(» N Sy\i£) = J N 2 Clf t \ogf t d^ +J N 2 C*J t d^ =J f t N 2 C N logf t dv» 



N 2 / f t (C N log f t 



£-N.ft 

ft 



)dv% + N z I C N f t dv 



.N 



Note that the last term is the price paid for not using an invariant measure. 

Since for every a, b > 0, a\og(b/a) — (b— a) is less than or equal to — (Vb— y/a) 2 , for every x, y € D N , we 
have 



ftC% v ,v log ft - c e * y Jt < -Pn{v ~ x, v) [v/M^o - VW) . 

An analogous calculation for the other parts of the generator permits to conclude that 



N 2 / f t (C N log f t 



C-Nft 
ft 



)du% <-N 2 D„.{f t ). 



To conclude the proposition we need a bound for N 2 j Cn ftdv^ ■ Let us write it explicitly: 



N' 



J £ N f t dv» = N 2 J (L^ft + C^ 2 f t + C c N f t + C b N f t )d^. 



Now, we compute each term inside this integral separately. 



TV 2 f C^hdv" = iV 2 jY, E Y[f(v-^, v +^ +ej , v )-fmd^ 

v£V x <zD d N j=l 

/d 
E E E[/^-^+^-^)-/w^ 

where d x . v represents a configuration with one particle at position x and velocity v, and no particles elsewhere. 
Then, if we let 

Jx , v = e v (A(K(x)))/(i-e v (A(K(x)))), 

the change of variables t] — V xv + d x + ej , v = £, changes the measure as dv^ (rj)/dv^ (£) = r ) x ^/l x +e j .v Hence, 
after changing the variables, we obtain 

d 



N- 



[c^f***? = ^ 2 EE / E 



+ ^EE / E 

veVj=l J x£D% 



= EE/ E 
+ ^E/E 



7:k,i> 

Aat7(x,i;) 



- 1 



- 1 



7s v 



JV 



si- 1 

^E/ E 

xi=iV-l 



1fx,v 



N 



Since 7 XjU is smooth and does not vanish, we can bound the above quantity by C\N d , where C\ is a 
constant depending only on a and /?. By a similar approach, one may conclude that 



Af/W<E& E 



which is clearly bounded by C^A^, where Ci is a constant depending only on a and j3. 

We now move to the generator with respect to collision. The change of variables r/ y,q = £ changes the mea- 
sure as dv™ {r])/dv£ (£) = (7^ 

,v'7y,w')i where v+u; — u'+u;'. Then, clearly, (jy, v Jy,w 

1, and therefore 



TV 2 1 C° N f t dv?=0. 



Lastly, we note that the change of variables a x ' v rj — £ changes the measure dv^ (r])/du^ (£) = a v (x/N)/(l— 
a v (x/N)) or (1 — a v (x/N))/a v (x/N), depending on whether there is or there is not a particle at the site x 
with velocity v, and analogously for (3. Therefore, a simple computation shows that 



N 2 J C h N f t dv»=Q. 



which concludes the Proposition. □ 

Let < /, g >„ be the inner product in L 2 (v) of / and g: 



< f,9 >v= J fgdv. 



Proposition 3.2. There exist constants C\ > and C2 = Ci(a, [3) > such that for every density f with 
respect to V^, then 

<C N ^f,^f > v n< -dD^(f) + C 2 N d - 2 . 

Proof: A simple computation permits to conclude that D^ N and D b N are both non-negative. Finally, the 
computation for Dffi follows the same lines as those on the proof of Proposition 13. 2\ and on Lemmas 13.41 
and!3.5i and is therefore omitted. □ 



3.2. Replacement lemma for the boundary. Fix k = 0, . . . , d, a continuous function G : [0, T]xT d 1 
" 1 , and consider the quantities 



Vk ( s > V> a i G) = jj-jZi Y G k(s,x/N)(l k (T) _^(s)) - ^2v k a v (x/N]j, 

V+(s,n,p,G) = J ^ rT G k (s,x/N)[l k ( V{N _ lti) (s))-Y / ^ v (x/N)), 

1 1 Ne ~ 1 

Vg(s,rj,a,G) = -^=1 Y G k (s,x)(j k (r} {hi) (s)) - — ^ 4 (?7(i,x) («))), 



reTl 



and 



X 1 = 1 



N-l 



Vk(s,r],P,G) = -^-j Y G k{s,x)(h{v(N~i, x )(s)) - ^j- Y h(ri(N-i,x){s)) 

xeTp 1 x 1= N(l-e)-l 

where s € [0, T], and G k , < k < d are the components of function G. 



The main result of this subsection is the following Lemma: 
Lemma 3.3. For each < t < T, < k < d, and G : [0, T] x D d 



continuous, 



lim sup E„ 



dsVl(s,r,X,G) 



= 0, 



where j — 1, 2, and £ = a, /3. 

Proof: It is clear that Vj! is bounded for each < k < d, and j = 1,2. By the entropy inequality, 



E„ 



d S Vi(s, V ,C,G) 



< 



HU N \ V %) 1 

AN d AN d B ^ 



exp 



dsAN d V£(s, V ,t,G) 



for all A > 0. We have that the first term on the right-hand side is bounded by CA , for some con- 
stant C. To prove this result we must show that the limit of the second term is less than or equal to as 
N — > 00 for some suitable choice of A > 0. Since e' x > < e x + e~ x and limsup^r^QQ N~ d \og{aN + on} < 
max{limsup Ar _^ 00 N~ d log(ajv), hmsupjv^^ N~~ d log(6jv)}, replacing Vj! by — Vji, or more precisely, replac- 
ing G k by — Gfc, we are able to conclude that we only need to prove the previous statement without the 
absolute values in the exponent. Let W k (s) = AN d Vj! (s, r/, G). Then, by Feynman-Kac's formula (see, for 
instance, [U03]), we nave 



E„ 



e^y\sAN d V^( Sl nX,G)^ 



where S^t is a semigroup associated to the operator = C + W k {t), for more details see [HI A. 1.7], see 
also [2]. Then, by Cauchy-Schwarz 

<5 ^i,i>^<<^n,^i>y N 2 . 



On the other hand, since W k is bounded, the adjoint in L 2 {v^) of Cf , Lf'* , is equal to C* N + W k {t). We 
have that 

= 2 < Cf* sZ k l, sZ k 1 >< A Wfc (s) < SZ" 1, S# 1 , 
where Avi/ fe ( s ) = su P||/|| i2( n } j { < ^fel 5 );/ >^ + < >i/™}- Therefore, we obtained that 

1 



log£^ 



exp • 



{| / d S AiV rf ^( s ,77,C,G)|} 



< 



N 



AN d ~ 2 



In this formula the supremum is taken over all densities / with respect to , and recall that < /, g >„ 
stands for the inner product in L (v) of / and g. An application of Proposition 13.21 permits to conclude 
that < CnVJi V7 >v n i s bounded above by CN d ~ 2 , where C > is some constant. Thus, if we choose, for 
instance, A = N, the proof follows from an application of the auxiliary Lemmas 13.41 and 13.51 given below. □ 



Lemma 3.4. For every 0<t<T,0<k<d, and every continuous G : [0, T] x T 



d-l 



lim sup E„? 

N— »oo 



daV£(8,ri,C,G) 



where £ = a, /3. 



Proof: We will only prove for a, since for the proof is entirely analogous. Note that G is continuous 
and its domain is compact, hence, we may prove the above result without G. Set f t = 1/t f Q f s ds. With 
this notation we can write the expectation above, without G, as 



h(V(i,£) ~ E V ka v {x/N) 



#5=r E f ftto) 

5=T E llM [v((l,5>),v)-a v (x/N)]d^. 



£6T d-i veV 



Then, splitting the integral into the integral over the sets [?7((1, x), v) — 0] and [77 ( ( 1 , x),v) — 1], and changing 
the variables as 1 — T)(xn, v) = £, we obtain 



N d - 



- E JfM 



h(v(i,x) - E v kOt v {x/N) 



dv, 



N 



where 



P Qi „ = a„(£/JV)(l - 77((1, 1), «)) + (1 - a„(i/iV))rK(l, 5), «). 

Writing {a — b} = {f t (r)) — f t (r) — Un t x),v)} as — Vb}{^/a + s/b} and applying Cauchy-Schwarz, the 
above expression is bounded by 

2t E ve v v k t - 
J + J^T^yffMth 

where D u iv b (f t ) is the Dirichlet form of f t with respect to C b N . Then, choosing A = y/~N, the proof of 



the Lemma follows from an application of Proposition 13.21 together with the fact that the Dirichlet form is 
convex. □ 



The next Lemma concludes the boundary behavior of the particle system. 
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Lemma 3.5. For each 0<t<T,0<k<d, and continuous G : [0, T] x D d , 



lim sup lim sup E t 

e^O JV-»oo 



N 



dsV k 2 (s,r),(,G) 



0, 



where £ = a, /3. 

Proof: First of all, note that since G is continuous and its domain [0, T] x D d is compact, it is enough to 
prove the result without the multiplying factor G. Moreover, we will only prove the first limit above, since 
the proof of the second one is entirely analogous. Considering the notation used to prove Lemma 13.41 we 
may write the expectation above, without G, as 

Ne—1 



t 



t E / J k(*i(i,s)) - E J * (*?(«!.«) 



dv 



N 



We now obtain, by a change of variables and a telescopic sum, that the absolute value of the above expression 
is bounded above by 



t 



Ne—1 y-1 



N d ~ 



y=l x 1 = l 



x\ — 1 



ft(U^M)-ft(l[ r z M) 



dv 



N 



where K\ is a constant which depends on a, (3 and d, z\ = 1, . . . , z y -\ = y is the path from the origin to y 
across the first coordinate of the space, and t Zi (77) • • • r Zi (77) is the sequence of nearest neighbor exchanges 
that represents the path along z±, . . . , Zi. By Cauchy-Schwarz, this expression is bounded above by 

2 



tA 



Ne-l y-1 



N d ~ 



t 



AN d - 



y=l X; — l 

Ne-l y-1 

Ne 



\ 



1=1 



\ 



XX— 1 



1 Ne-l y-l 



xEV' 



y=i x 1= l 



1=1 

xi—l 



dv 



N 



ft([[T z M)-f t (U r z M) 



i=l 



i=l 



dv 



N 



for every A > 0. Now, we can bound above the last expression by 

tAKl TV* (1 \ -i- ^ aiVe 
J^i^u^Ut) H ^ - 

for every ^4 > 0, where if 2 is a constant that depends on K\. Then, choosing A = y/iN and applying 
Proposition 13. 21 we conclude the proof of this Lemma. □ 



3.3. Tightness. To prove tightness of the sequence (Qn)n, it is enough to prove that for every k = 0, . . . ,d 

= 0, 



lim limsupE^N 

<5^0 n^oc 



sup 

\t-s\<8 



i& E H (^) '*(*(*))- ii E 

xeD% x£D d N 

for any smooth test function H : D d — > R vanishing at the boundary. 
Fix < k < d, then, by Dynkin's formula 

11 irk ^ k,N u . . fc.JV _ f r k N tt j 

M t =< 7r t ' , H > - < ir a > — £n < ,H > ds 

Jo 

is a martingale. On the other hand, 



/ {£at < tt^ n ,H> 2 -2 < n k s ' N ,H> C N < n*< N ,H>}di 
Jo 



(3.3) 



Writing the above expression as four sums, the first corresponds to the nearest neighbor symmetric exclusion 
process and the other corresponds to the asymmetric exclusion process, the third and fourth corresponding 
to the collision and boundary parts of the dynamics, respectively. A long, albeit simple computation shows 



that all of these sums are of order (D(N~ d ), and therefore, the right-hand side of the above expression is of 
the same order. Thus, by Doob's inequality, E^n [sup 0<s<t (M t fe ) 2 ] = 0(N~ d ). 
Hence, by (|3.3p and the above estimates, we have 



w E 5 (J)w)) = f E h(^) + 

1 ^ /"* x 

+ JjdYl E El p(z,w)wfery r (0,w)[l - ry r (z,w)]z 3 (a Uj iJ) + 

pE /W)(|)4fe(r))dr + ^ T 2 f d Ul H I k (rj x (r))dr 

El=JV-l 

£ ^ t a illJ ff(^)/ fc fe(r))dr + J R Ar + 0(iV-' i ) + 0(iV- 1 ), 



xeD% 

xi—l 

where the terms were obtained from £n < Ttg' N , H > by means of summation by parts, and the replacement 
of discrete derivatives and discrete Laplacian by the continuous ones, and Rn is the error coming from such 
replacements. Since p is of finite range, the error Rn is uniformly of order 0(N^ 1 ). Finally, by using Lemma 
13.31 and a calculation similar to the one found in equation (|3.10p . we have that C b N < Kg , H >= 0(N^ 1 ). 
Tightness thus follows from the above estimates. 

Our next goal is to prove the replacement lemma. To do so, we need the following result known as 
equivalence of ensembles, which will be used in the proofs of the one block estimate and of the two block 
estimate. 

3.4. Equivalence of ensembles. FixL> 1 and a configuration rj, let I L (x, rf) :—I L {x) — (Iq (x) , . . . , 1% (x)) 
be the average of the conserved quantities in a cube of the length L centered at x: 

1 L < zex+A L 

where, Al = {—L, . . . , L} d and |Al| = (2L + l) d is the discrete volume of box Al- 

Let 2Jl be the set of all possible values of I L (0, rj) when r/ runs over ({0,1} V ) , that is, 

Z3 L = {l L (0,r 1 y,r,e ({0,1} v ) Al }. 

Note that QJl is a finite subset of the convex envelope of {/(£) : £ G {0, 1} V }- The set of configurations 
({0, 1} v ) Al splits into invariant subsets: for i in 23l, let 

n L (i) := {ve ({0,l} V ) Ai :I L (0) = i}. 

For each i in 03^, define the canonical measure i>L,i as the uniform probability measure on TL^i). Note that 
for every A in M d+1 

"A*,iO = (•| /L=i )- 

Let < g; f >fi stands for the covariance of g and / with respect to /i: < g; f > M = E^[fg] — E^[f]E^[g\. 

Proposition 3.6. (Equivalence of ensembles): Fix a cube Kg C A^. For each i E %3l, denote by ir the 
projection of the canonical measure VA L ,i on A^ and by \i the projection of the grand canonical measure 
A'A(i) on Then, there exists a finite constant C(£, V), depending only on I and V , such that 

\EAf]-EAf]\<^<f-J>li 2 

for every function f : ({0, l} v ) Af i— > K. 

The proof of Proposition 13.61 can be found in Bcltran and Landim [lj. 
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3.5. Replacement lemma. We now state the replacement lemma that will allow us to prove that the limit 
points Q are concentrated on weak solutions of (|2.4[) . 



Lemma 3.7. (Replacement lemma): For all 8 > 0, 1 < j < d, < k < al: 



lim sup lim sup I 

e^O N^,oo 



[ £ T x Vg( V (8))d8>5 



= 0, 



where 

vt k (v) = 



(21 



y&A e v£V zi 



VEV 



(3.4) 



Note that V^j? is well-defined for large N since p(-, v) is of finite range. We now observe that Propositions 



and 13.21 permit us to prove the following replacement lemma for the boundary driven exclusion process 
by using the process without the boundary part of the generator (see [17] for further details) . We postpone 
the rest of the proof to Section |U 

3.6. Energy estimates. We will now define some quantities to prove that each component of the solution 
vector belongs, in fact, to iJ 1 ([0,T] x D d ). The proof is similar to the one found in [15] . 
Let the energy Q : D([0,T],M) — > [0, oo] be given by 

d 

Q( 7r )=^Q i ( 7r ) J 



with 



sup 

Gec~(n T ) I Jo 



2 / dt < n t ,d Uz G t > 



duG(t, u) 



o J D d 



where f2y = (0,T) x D d and C£°(f2:r) stands for the set of infinitely differentiable functions (with respect 
to both the time and space) with compact support in fi^. Let now, for any G S C^°(CIt), 1 < i < d and 
C > 0, Qf c : D([0, T],M) — > K be the functional given by 

Q?c( 7r )=[ ds <n s ,d Ul G s > -C [ ds [ duG(s,uf. 

JO Jo J D d 



Note that 



sup {Qf tC } 
Gea~(n T ) 



4C 



Lemma 3.8. There exists a constant Co = Co (ft) > 0, such that for every i = 1, . . . , d, every k = 
and every function G in C£°(fir) 



(3.5) 
d. 



) " ' ) ■*) 



Um sup _L bg E v n [exp{iV d Qf Co (^ fc )}] < C . 



log E v * 



expJiV^ ds {I k {r] x (s))-I k {r] x - ei {s)))G{s,x/N) 



Proof: Applying Feynman-Kac's formula and using the same arguments in the proof of Lemma 13.31 we 
have that 

j_ r r ,t 

N d 

is bounded above by 

where is equal to 
sup 



1 

N 1 



\?ds, 



ap{(rt%I k (ri(x)) - h(v(x - ei)))G(s, x/N)J^ + N 2 <£ N ^f,yff >„„ }, 
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where the supremum is taken over all densities / with respect to v£ . By Proposition 13.21 the expression 
inside brackets is bounded above by 

N 2 

T 



CN d -—D v s(f)+ \NG(s,x/N) I [I k (Vx) ~ HVx-eiMvKidv) 

XED* 



We now rewrite the term inside the brackets as 

E Ufc E { [ NG( S ,x/N)[r 1 (x,v)- V (x-e l: v)}f(r 1 )^(df 1 )}. 

After a simple computation, we may rewrite the terms inside the brackets of the above expression as 

NG(s, x/N) [ [r,(x, v) - V (x - e h v)]f(r,)^(d V ) 



(3.6) 



= NG(s,x/N) J V (x,v)f( V )^(dr 1 ) 

NG(s,x/N) f v (x,v)f(r ] x - e ^) 1 ^ i ^(dv) 

J Ix.v 

= NG(s,x/N) I v(x,v){f( V )-f( V *~ e ^n}^(dv) 

'fx — Ci ,v 



+ G I r](x,v)f(r] x ~ ei ' x ' v )N 



1 



< G{s.x/Nf I f( V *- e ^)^(dv) 



IX, V 



+ 



r}(x,v)f(rf -«>*•«) 



Nil- 



1x — ei i V 
^yx,v 



+ 2G(s, x/N) 2 J V (x, u)(v7fo) + V ' f{rf-*>*>«)) 2 V » (d V ). 

Using the above estimate, we have that (|3.6|) is clearly bounded above by C\ + CiG(s,x/N) 2 , by some 
positive constant C\ = C\{k), since j. tV is smooth and the fact that / is a density with respect to . Thus, 
letting Co = C + C\, the statement of the Lemma holds. □ 

It is well-known that Q(tt) is finite if and only if it has a generalized gradient, V7r = (d Ul n, . . . ,<9 Ud 7r), 
and 

Q(tt) = f [ du\\X7Tr t (u)\\ 2 < oo. 



J D d 

In which case, Q(tt) = Q(tt). Recall that the sequence (Qn)n defined in the beginning of this section is 
tight. We have then the following proposition: 

Proposition 3.9. Let Q* be any limit point of the sequence of measures (Q n )n- Then, 



and 



E Q . 



En* 



ds / \\Vp(s,u)\\ 2 du 

n \JD d 



ds / ||Vpfc(*,u)|rdu 

\JD d 



< OO, 



< oo. 
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Proof: We thus have to prove that the energy Q(n) is almost surely finite. Fix a constant Co > satisfying 
the statement of Lemma [3.81 Let {G m : 1 < m < r} be a sequence of functions in Cq^SIt) (the space of 
infinitely differentiable functions vanishing at the boundary) and 1 < i < d, and < k < d, be integers. By 
the entropy inequality, there is a constant C > such that 



E„ 



max 

Km<r 



<C+±- d \o g E v 



exp 



^nmx fe^)} 

Krn<r I ' u ) 



Therefore, Lemma 13.81 together with the elementary inequalities 



limsupiV d log(ajv + 6jv) < lim sup maxjlim sup N d log(aAr),limsup N d log(&Ar)} 



and exp{max{ 



Xi, 



i}} < exp(xi) 



exp(x„) imply that 



E Q . 



max 

Km<r 



{e&(***)} 



= lim E.. 



max 

1< m < r 



Kc (^' fc )} 



Using this, the equation (|3.5p and the monotone convergence theorem, we obtain the desired result. □ 



3.7. Proof of Theorem 12.11 Note that all limit points Q* of (Qn)n are concentrated on absolutely 
continuous measures with respect to the Lebesgue measure since there is at most one particle per site, that 
is, 

<3*{tt; ir k (du) = p k {u)du, for all < k < d} = 1, 
where n k denotes the fcth component of tt and p = p. 

For k — 0, . . . , d, denote, again, by 7rf ' N the empirical measure associated to the fcth thermodynamic 
quantity: 



Jx/N- 



k,N,b x 



and 7r. 



Further, denote by n t 
quantity restricted to the boundaries: 



k,N,b N - 



1 the empirical measures associated to the A:th thermodynamic 



k,N,bi 



N d 



l^T h{vx(t))S x/N , 



x£D d N 



for i = l,N-l. 

To compute £n < 7r t ,JV , H > for this process, we note that £ L N Ik(Vx) vanishes for k = 0, . . . , d, because 
the collision operator preserves local mass and momentum. 

Since, in our definition of weak solution we considered test functions H vanishing at the boundary, that 
is, H{x) = 0, if x € {0, 1} x T d_1 , we assume that H vanishes at the boundary as well. 

Now, we consider the martingale 



M, 



N,H 
t.k 



~-<n t ' ,H>-<n ,H> 



N 2 C N < ir*> N ,H> ds, 



which can be decomposed into 



N,H 
t.k 



, k,N tj . . k.N u . 

<TT t ,H> — < 7T ,H> 



2 rex,l 



N Z C 



< n*< N ,H> ds 



1 r ex,2 



< tt^ n ,H >ds- I N' Z C° N < ir«'" ,H>ds. 



N 



2 r b 



k,N 



We first prove that 



N 2 £% < TT 



k,N 



H > ds 



(3.7) 
(3.8) 

(3.9) 



vanishes as N — > oo. A simple calculation shows that 

N 2 £ b N r)(x, v) = N 2 [a v (x/N) - v {x, v)] , if xi = 1, 
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and 

N 2 C b N 7](x, v) = N 2 \j3 v (x/N) - r)(x, v)} , if x± = N — 1. 

Since H vanishes on the boundary, H((x + e\)/N) — if x\ = N — 1, and H((x — ei)/N) = if x\ =0. 
Then, we have the equalities NH(x/N) = d*H((x - e^/N), if Xi = 1, and NH{x/N) = -d*H(x/N), if 
x\ = N — 1. Therefore, we obtain 

N 2 C b N <^ N .H> = 3 ^=tEx€^E 06 v«*[«»(#)-'?(^«)]^^(t^) 

- l^ET e W E, e v^[/?(#)-^«^(f). (3 ' 10) 

We now use the last computation together with Lemma f3T3l to conclude that (|3.9j) vanishes as N — > oo. 
Further, after two summations by parts of the integrand on the right-hand term of (|3.7|) . we have that 

N 2 C e ^ <^ N ,H >ds = I f <n*< N ,A N H >ds 

* Jo 



and after one summation by parts on the right-hand term of (|3.8p . and noting again that H vanishes at the 
boundaries, we have that 



t ^ i-t d 



N 2 £ e *' z <ir k s > N ,H>ds 



J0 i=ixexi 



where t x stands for the translation by x on the state space Xm so that (T x ri)(y,v) = rj(x + y,v) for all 
x, y € Z d , v e V, and W^- is given by: 

= ^2 vk Yl p(*> u )zj^(°> u )[i-»7s(*,i>)] ) 



where wq = 1. Since v) is of finite range, 



where x( a ) — a (l — a )- Now, note that E v n (r)(x, v)) = a v (x/N) if x G {1} x 1 and E u n(t](x,v)) 
Pv{x/N) if x € {N - 1} x T^ 1 . 

We then apply Lemma 13.71 to write the martingale in terms of the empirical measure. Further, we apply 
the replacement lemma for the boundary (Lemma 13. 3j) to obtain that all limit points satisfy the integral 
identity in the definition of weak solution of the problem (|2.4p . 

Using the previous computations and the tightness of the sequence of measures Q n (for more details see 
[T4l Chapter 5]) we conclude that all limit points are concentrated on weak solutions of 

d t (p,p) + « [« • v x (e v (A(p,p)))\ = \a(p, P ), 

with boundary conditions, given in the trace sense, by 

(p,p)(t,x) = a(x), for x £ {0} x T d_1 , (3.11) 

and 

(p,p){t,x) = b(x), for x G {1} x T d ~\ (3.12) 

where a(-) and &(•) were defined in equation (12. 3j) . and vq — 1. The uniqueness of weak solutions of the above 
equation implies that there is at most one limit point. Moreover, by Proposition 13.91 each limit point of 
(Qn)n is concentrated on a vector of measures with finite energy, that is: whose components have densities 
with respect to the Lebesgue measure that belong to the Sobolev space H 1 (D d ). This completes the proof 
of the theorem. □ 
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4. Proof of the replacement lemma 

As mentioned in the subsection [331 we only have to prove this result for the process without the boundary 
dynamics. In this case, we have a product invariant measure given by v^ p . 
Let Ijl n {T) be the Cesaro mean of fi N S^ , namely: 

and let f Tk be the Radon-Nikodym density of fi N (T) with respect to . We have that the Dirichlet 

form of f T ^, DN(fT.ki v p. P )i is bounded by CN d ~ 2 /2T, where C is some constant. Therefore, to prove the 
replacement lemma, it is enough to show that 

limsuplimsup sup / —: V T x V${r){s))f(ri)v" (dq) = Q. 

e ^ N^oo D N (f,u p , p )<CN"-^ J ^ x ^ 

From now on we will simply write the Dirichlet form of a function / with respect to the measure v^ p as 
D N (f). 

To prove the replacement lemma, we will prove the one and two block estimates: 
Lemma 4.1. ( One block estimate): For every constant C > 0, for 1 < j < d and for < k < d: 

limsuplimsup sup / V (T x V e 3 ' k )(i])f(i])Vp p (dri) = 0, 

t^oo JV-too D N { f)<CN d ~ 2 J iV 



where V^' k (ri) was defined in Lemma \3. 7[ 

Lemma 4.2. (Two block estimate): For every constant C > 0, for 1 < j < d and for < k < d: 
lim sup lim sup lim sup sup sup / — —? \l e (x + y) — I Ne (x) f(r])i/^ = 0, 

£-»oo JV^oc D N (f)<CN d - 2 j/SA.n J N ~[ d 1 



where I (x) was defined in subsection^, 

4.1. Proof of one block estimate. We begin by noting that the exclusion rule and the fact that V is finite 
prevents large densities or large momentum on I e (0). 

We have that the measure v^ p is translation invariant. Therefore, we can write the sum on one block 
estimate as 



V, 



N d 



where / stands for the space average of all translations of /: 

Denote by X( the configuration space ({0, 1} V ) Ae , by £ some configuration on X( and by iA the product 
measure v^ p restricted to Xf,. For a density / : X^ — > M + , fa stands for the conditional expectation of 
/ with respect to the er-algebra generated by {rj(x, v) : x £ A(,v £ V}, that is obtained by integrating all 
coordinates outside this hypercube: 

fl{Xi) = v ( 77% J 1 {r r .r,(z,v)=£(z,v), Z <£A i ,v£V}f( r l)Vp,p(dri), 

for e € X e . 

Since V/' k (r]) depends on the configuration r\ only through the occupation variables {t](x,v) : x £ A.£,v £ 
V}, in the last integral we can replace / by f e . In particular, to prove the lemma it is enough to show that 

limsuplimsup sup f V?'*(07*(0»£p(d0 = °- t 4 - 1 ) 

>oo JV^oo D N (f)<CN d - 2 J 
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We will now compute some estimates on the Dirichlet form. Let < •, • > v be the inner product in L 2 (v). 
For positive /, denote the Dirichlet form of f as: 

D N (f) = -< x/f,(£%+£ c N )f> <%> 

:= D NA (f)+D N . 2 (f)+D N .Af)- 



We have that 



and 



where 



d n M)= E j #(/)> 

x.zeD% 

\x — z\ — l 
x,zeD% 

D N M) = E 4 c) (/), 



v£D d N 



4%f) = E \ [ wsijf M ) y/m\ 2 < P (dv), 

4 2 l(f) = E / p(z,v)Wf(v xM )- x/W)M P ( d v) 



and 



4?K.f) = E / - mMpW- 

qeQ J 

Since the Dirichlet form is translation invariant and convex, we have that -Dat(/) < Djv(/)- 
Now, let 

E W) + E ^ a) w+E J - ,(c) w. 

|x-z|=l 

where each equals /W with i/^ p replacing z/^ p . By using Schwarz inequality and the definition of fi, 
we obtain that 

iifHJe) < 4]l(7),iifH7e) < 4 2 1(7) and i^) {lt) < j(c) (7) 

for every i,zeA(. Therefore, 

^U)< E 4:1(7.) + E ^4 2 1(7,)+E4 C) (L). 

x,z£Ai x,z£Ai x£Ae 

\x-z\=l 

On the other hand, by translation invariance of/, I x ]l(f) = Ix+y tZ+y {J), IxJiJ) = if+y^+yQ) and I x c) (f) = 
4 C) (7). Hence, 



D'Qi) < (2t + ir e #i (7) + E 42 (7) + (a + i) d 4 c) (7) 

i 

(2£+l) d 



N 

yeA,. 



< 



^(^.ll/) + Av, 2 (/) + Av,c(/)). 



Since the Dirichlet form is positive, £>jv(/) < CW d - 2 implies that D NA {f) < CN d - 2 , D Na {f) < CN^ 1 
and D N>c (f) < CN d - 2 . Thus, 

D e (fi) < 3C(2£+l) d N- 2 := C (C, ^iV" 2 . 
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Therefore, the Dirichlet form of f e vanishes as N f oo. Hence, by (|4.ip . to prove the one block estimate we 
must show that 



limsuplimsup sup [ V?" (QfiQv 1 (d£) = (4.2) 

£—00 TV— 00 D e (f)<C (C,i)N- 2 J 

with the supremum carried over all densities with respect to v 

We wi 
therefore 



£—00 TV— 00 D l (f)<C (C,l)N- 

,N 

We will now take the limit as N t 00. To do so, we note that V/' fc < C\, where C\ is some constant, and 



This subset of A4+(Xe) is compact for the weak topology, and since it is compact, for each TV, there exists a 
density /jy with Dirichlet form bounded by CqN~ 2 that reaches the supremum. Let now N n be a subsequence 
such that 

lim f y/'Viv„(0< P (^)=limsup f V^OM^p, P (^)- 

To keep notation simple, assume, without loss of generality, that the sequences N n and N coincide. By 
compactness, we can find a convergent subsequence /iv n . Denote by /oo the weak limit. Since the Dirichlet 
form is lower semicontinuous 

£>*(/oo) = 0. 

Moreover, by weak continuity, 

In conclusion, expression (14. 2[) is bounded above by 

lim sup sup / Vt\0f{04, P {di). 

l^oo D e (f)=0J 

We will now decompose along sets with a fixed number of conserved quantities. 

Recall that 53 l is the set of all possible values of I (0) when r\ runs over ({0, l} v ) Ai . Further, 5Jl is 
finite. Furthermore, consider for each i in QJl the canonical measure VL,i defined in subsection 13. 4| and 
moreover, recall that 

"jwO = /4 1 (•| jL = i )- 

A probability density with Dirichlet form equal to zero is constant on each set with a fixed number of 
conserved quantities. It is convenient therefore to decompose each density / along these sets. Thus 

V/'*(0/(0<p(dO = E W) I V^vlM), 

where, 

Since J2jeH e (j) Tj(f) = 1, to conclude the proof of the one block estimate, we must show that 



limsup sup / Vf' (£)v£j(d£) = 0. 
£^00 jezii J 



Since the measure vij is concentrated on configurations with conserved quantity j, the last integral equals 



/ foe 1 T)d E J2 Vk J2 p ( z > v } z i T y( h (^ z ' v ))-^ v J v k E AHteuvfi 



where h(£,z,v) = £(0,v)(l - £(z,v)). 

Fix some positive integer n, that shall increase to infinity after I. Decompose the set Kg in cubes of length 
2k + 1. Consider the set A = {(2n + l)x, x S Z d } n A^_„ and enumerate its elements: A = {xi, . . . , x q } in 
such a way that \xi\ < \xj\ for i < j. For 1 < i < q, let i?.; = + A„. Note that fl Bj = if i ^ j and 
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that Ui<i< q Bi C Ag. Let Bo = Ag — Ui<i< q Bi. By construction \Bq\ < Kn£ d 1 for some universal constant 
K. The previous integral is bounded above by 



veV \ l ' y£Bi z ' J 



Since \B \ < Knl d 1 , Y, v Ufc£(0, v)(l-£(z, v)) has mean Y, v «fcX(^«(A(j))), and l^eBi^C ' is bounded, 
the sum is equal to 



E^ [jbi E Z>( 



z, v)z j T y {h{£ > , z, v)) - VjE u t [h(Z, et,v)] 



plus a term of order 0{n/tj. Since the distribution of {£(z,t>); z 6 £ V} does not depend on i, the 

previous sum is equal to 



/ E Wfc I (o +i)d E '52p( z > v ) z 3 T y( h (t' z ' v )) - l 'j E A h (tei,v)} 



plus a term of order 0(nj€). 

Now, let fi\ be the product measure on ({0, 1} V ) with marginals given by 

fi x {v ■ r)(x, •)=£} = m A (£), 

for each £ £ {0, 1} V and at G Z d . Therefore, £L/[f(0,u)(l - £(euv))] = [£(0, u)(l - £(ei,u))l, where 
Vj = HAtj)- Moreover, if in the equivalence of ensembles we choose L — L(£) — [C(£, V)J, where C(£, V) 
is the constant given in the equivalence of ensembles, we can replace the canonical measure by the grand 
canonical measure paying a price of order og(l). Therefore, we can write the previous integral as 



/ E W M To ZTw E *52p( z > v ) z i T y( h (^ z , v )) ~ VjE Vj [h(£,ei,v)] 



plus a term of order og(l). We now note that Vj equals ir. on Ag . Then, the integral can be written as 

it now, 

V ^eA„ z / 



plus a term of order o^(l). Let now, 
but we know that ei , t>)] = x(0„(A(,j))), then 



vev 



9i (0 



E 



(2n 



Now, ({0, 1} V ) is compact on the product topology, and also, all the marginals of Vj converge to the 
marginals of v p , p , when j — > (p,p) as I — > oo. Then, converges weakly to ^ p . p . Further, since gj(£,) — > 
5p,p(£) for every £, we have from Theorem 5.5 of Billingsley [6], that 



£->oo 
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this convergence being uniform on compact subsets of R+ x M. d . Then, since the remainder term is 01 (1), 
the limit as £ — > oo and j — > (p, p) is 



/ (2n+ l) d ^ ^ v k^ z jP{ z i v ) T y(K£,i z iV))-^v 3 v k x{0v{k{PiP))) 



v p ,p{d£). 



On the other hand, as k 1 oo, by the law of large numbers, this integral converges to 0. 
Therefore, the one block estimate is proved. □ 



4.2. Proof of the two block estimate. To prove the two block estimate, it is enough to show that 
lim sup lim sup lim sup sup sup f — ^ E I^O") 

X ^ N 

-l\x + y)\f( V )v» p (di 1 )=0. 
As for the one block estimate, we can rewrite this integral as 

I<(0) -I\y) 7(vK P (d V ), 



(4.3) 



/ 



where / stands for the average of all space translations of /. I e (0) and I e (y) depend of the configuration i] 
only through the occupation variables {rj(x, v) : x G A Vt £, v S V}, where 



A,, 



{-£,..., £} d U[y + {-£,...,£}" 



We now introduce some notation. For positive integer £, let X 2 - 1 denote the configuration space ({0, 1} 
({0, l} v ) Af , £ = (^1,^2) the configurations of X 2 ' e and the product measure v^ p restricted to X 2 ' 1 (which 
does not depend on N) by v 1 '„. Denote by f V} e the conditional expectation of / with respect to the tr-algebra 
generated by {rj(x,v) : x <G Ay,i,v 6 V}. 

Since I e (0) and I^(y) depend on i](x, v), for x G A Vt £ and v G V, we may replace / by /„ e , and then, we 
can rewrite f|4.3[) as 

lim sup lim sup lim sup sup sup / V l-E^(O) - E e 2 (0) f y A€)uZ'*(d£) = 0, 

£->oo e^0 JV^oo D N (f)<CN d - 2 y£(A eN \A e ) J ^ TT., 1 



where 



E 



1 1 zG^+Ai 



Now, we need to obtain information concerning the density f yt from the bound on the Dirichlet form of /. 
Then, let D 2 ' e be the Dirichlet form defined on positive densities h : X 2,i — > E + by 



where, 



^w=E / 

v£V J 

+ EE/ 



, 2 TV 

\x-z\=l 



vGV J 



E ? + ^7 E 



, 2 TV 

|x-z|=l 



and, 



E 9+ 



+ E 



_i*i=i 

E / p(°>9,£i) 



H£'~' v ,£' + ' v ) - VW) 



i>GV 



1 1 / X 



1*1=1 



where 



£i(0, v) ±1, if a; = and w = v, 
£i(x,w), otherwise. 
This Dirichlct form corresponds to an interacting particle system on (V x Af) x (V x A|), where particles 
evolve according to an exclusion process with collisions among velocities on each coordinate and where 
particles from the origin of one of the coordinates at some velocity can jump to the origin of the other at 
this velocity and vice- versa. 

Using the same idea as for the one-block estimate, we can prove that 

D{(J y j) < D N (J) and D e 2 (f y/ ) < D N (f), 

and hence, 

D[(f y , t ) + D l 2 Cf yA )<2CoN-\ 
for every density / with Dirichlet form Djq{f) bounded by CN d ~ 2 . It remains to be shown that we can also 
estimate the Dirichlet form io,o(/j/,£) by the Dirichlet form of /. 
We begin by noting that 

where, 



4:Uh) = E 



E + N P ^ V ^ 



1*1 = 1 



+ 



and 



J o,oW = E /MO, 9, 6) 
vev J 

+ E / P(o,<z,6) 



1 2 



Then, a simple calculation shows that 



4; 2 (/^)<24 c) (/), 



and therefore Iox>(f y ,e) is a l so 01 order N 2 . We then have to obtain a bound for Io o(f y e)- 

Following the same lines used to prove that ii'jPfft) < Mif) in the proof of the one block estimate, for 
j = 1, 2, c, we have that each density /, with respect to v^ p , /p'j(/ 5 ^), is bounded above by: 



1*1=1 



(4.4) 



Let (a;fe)o<fc<|||2/||| be a path from the origin to y, that is, a sequence of sites such that the first one is the 
origin, the last one is y and the distance between two consecutive sites is equal to 1: 

xq — 0, %\\\y\\\ = y and \xk+i — Xk\ = 1 for every < k < \\\y\\\ — 1, 

1 1 • 1 1 1 is the sum norm: 

|||(» ls ...,%)|||= E \Vi\- 
i<i<d 

Let r Xl • • ■T Xi (rj) be the sequence of nearest neighbor exchanges that represents the path along xi, . . . ,x%. 
Then, by using the telescopic sum 



j \ J ±1 

\ \ i=l 



/(II ^.M) 



\ 



k-1 



and the Cauchy-Schwarz inequality 



'111*111-1 



£ fulfill E 



we obtain that (|4.4[) is bounded by 
2£ 



E J + 



1*1=1 



insni-i 
IMII E 

fc=0 



fe=0 



\ 



/(ii T -.(i))- 



\ 



k-l 



/(II T *M) 



IIMIM 



fe=0 



Since / is translation invariant, for each k, l2h,x k +i if) = ^x^+z x k+l +z{f) f° r au z G Hence, /|2x fe+ i (/) < 
N^ d D^(f). In particular, 

4;o(7^)<2 d+2 |||y||| 2 ^- d J Div(/). 

Recall that ?/ G A e jy, and hence |y| < 2iVe, | • | is the max norm. Then, |||y||| < d\y\ < IdNe. Since the 
Dirichlet form is assumed to be bounded by CN d ~ 2 , we have proved that 

4;o(7^) < 2 d +U 2 Ce 2 . 

We have, therefore, proved that for every density / with Dirichlet form bounded by CN d ~ 2 and for every 
e?-dimcnsional integer with max norm between It and INe, 

D 2 ^(J y/ )<C 2 (C,d,i)e 2 . 

We can now restrict ourselves to densities / such that D 2 (f y g) < C^e 2 , that vanishes as e J, 0. In particular, 
to conclude the proof, it is enough to show that 

hmsuphmsup sup / \E[{Q) - E{{Q)\f ^)vW{d£) = 0, 

l^oo e^O D 2 - e (f)<C2(' 2 J 

this time, however, the supremum is taken over all densities with respect to v 2 ^. The rest of the proof 
follows the same lines as the ones in the one block estimate, beginning by decomposing the Dirichlet form 
along the sets having fixed conserved quantities and then applying the equivalence of ensembles. Therefore, 
the two block estimate is proved. □ 
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5. Uniqueness 



To conclude the proof of the hydrodynamic limit, it remains to be proven the uniquenesses for the solutions 
of problems (|2.4[) and (|2.4p . The strategy we used to prove this result was employed by Oleinik and Kruzhkov 
[TO] and is due to Yu.A. Dubinskii. 

Let v and to be two weak solutions to the problem (|2.4[) . corresponding to the same initial function Uq. 
Fix some j = 1, . . . , d + 1, and let Hj € C 1 * 2 ([0, T] x D d ) be such that Hj (T, u) = 0, for all w. Then the 
integral identity for v — u> holds: 



dt 



D rl 



du(Vj — OJj) 



Ki<d 



h / eft / v i(9v(v) - g v {uj)) E v idu z Hj = 0, (5.1) 

• /o ^ «ev i<i<d 



where g v {v) = x(#„(A(j/))), v j^j an d Hj are the components of v,u> and respectively. If Vj = ojj, we 
already have what we want, thus, suppose Vj ^ u>j. Introducing the notation 



Pi 



3 _ 9v{v) -9v{u) 



we have that we can write (|5.1|) as 

r-T 



dt 



''"(.",■--•,■) w + 5 E ^ + E^ E v ^ H o =°- 8 5 21 

l<i<d veV l<i<d 

Now, let /3j m be a sequence of smooth functions which converge in L 2 ([0,T] x D d ) to ftl, as m — > oo. We 
denote by HJ™ (t, u) the classical solution of the equation 

d * H T + \ E d l H T + E «i^ ,Tn E V ^ H T = ^> ( 5 - 3 ) 

l<i<d vGV l<i<d 

H?(T,u)=0,H?(0,u) = 0, 

where $j is a smooth function finite in [0, T] x Z? d . For more details on the solutions of partial differential 
equations of the parabolic type, the reader is referred to Friedman [13] , and for details on solutions of systems 
of linear partial differential equations of the parabolic type in general, the reader is referred to Ladyzenskaja 
et al. [15] . 

Now, if we replace Hj in ()5.2|1 by HJ 1 and use (|5.3|) . we obtain: 



dt 



D'l 



du{vj -Uj)$j + 



dt 



vev 



Ki<d 



= 0. 



(5.4) 



Finally, since we are in a compact domain and the coefficients /^ ,m are smooth, we have that there exists 
an M > such that \H™\ < M. Since these coeffiecients converge in L 2 ([0, T] x D d ) 1 the constant M may be 
taken to be independent of m. Multiplying (|5.3p by HJ 1 , integrating over [0, T] x D d , and then integrating 
by parts, we have that 



fdtj du^U^f-) 2 = fdt f du lj2^Pi' mH T E v ^ h t - ® h t) - \f mh 

Jo Jd* i=1 * V ou, j j j Dd y vev ^ d j 2 j Dd 

On applying the elementary inequality \ab\ < ea 2 + 6 2 /(4e) and using that \H" l \ < M, we obtain that 



f dt f duj2l 

Jo Jd" £~T 1 



1 (dHJ 



V dui 



<C, 



where C is a constant that may depend on M and $, but not on m. 

Therefore, by applying the Cauchy- Schwartz inequality and using that j3i' m converges to (3 3 v in the L 2 - 
norm, we see that the second term on the left-hand side of equation (|5.4|) tends to zero as m tends to infinity. 
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This implies that for every e > there exists m such that the absolute value of the second term on the 
left-hand side of equation (|5.4[) is less than e. We, then, have obtained that 

Ve > : / dt du(vj - oJj)^j < e, 
Jo JT d 

and hence, for each j — 1, . . . , d + 1, Uj = u>j. Therefore v = u). □ 
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